A lattice study of the quark propagator and vertex function 
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We report on the status of a study of the quark propagator and quark-gluon vertex in momentum space. 
Quark propagators have been generated at (3 = 6.0 using the 0(a)-improved Sheikholeslami-Wohlert (SW) action 
and fixed to the Landau gauge. The first results for the quark pole mass and field renormalisation constant are 
reported, and plans for future work are presented. 
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1. Introduction 

The study of quark and gluon correlation func- 
tions on the lattice provides a powerful tool for 
increasing our understanding of QCD. Firstly, it 
enables us to study non-perturbatively the funda- 
mental degrees of freedom in QCD, ie the quark 
and gluon fields. Secondly, it allows for a more di- 
rect study of phenomena like confinement, hadro- 
nisation and dynamical chiral symmetry break- 
ing. Thirdly, it provides an opportunity to ex- 
tract basic QCD parameters like the running 
coupling and quark masses from first principles. 
Here, we will report on work that is in progress 
to study the quark propagator and quark-gluon 
vertex. 

2. The quark propagator 

The quark propagator may be studied both in 
time and in momentum space. The first approach 
is suited to investigating certain aspects of dy- 
namical mass generation, while the second ap- 
proach allows us to extract field renormalisation 
and the quark pole mass. 

2.1. Quark propagator in time 

The zero-momentum quark propagator in a co- 
variant (or related) gauge is expected to show the 
following behaviour [1]: 



s(t) = ($>0MMo,o) 

\ X I 

= A(t)+ 74 B(t) 



where A(-t) = -A(t) and B(-t) = B(t), with 
the asymptotic behaviour 



A(t) 
B(t) 



Ae ~m s t 
Be -m v t 



(2) 
(3) 



The 'scalar' and 'vector' quark masses m s and 
m v should be equal, as should the amplitudes A 
and B. 

2.2. Propagator in momentum space 

At tree level, the continuum momentum space 
propagator is 



i f + m 



(4) 



On the lattice, we would therefore expect the 
inverse propagator for some appropriate range of 
momenta (1/L <C p <€i l/ a ) t° behave according 
to 



r(p) = (pa) (i f + m(pa)) 



(5) 



where m(pa) is the pole mass, which can be com- 
pared with calculations in perturbation theory 
and from Schwinger-Dyson equations. 

3. Vertex function 

We can define a quark-gluon vertex on the lat- 
tice with quark momenta p — q and p, and gluon 
momentum q as 



(S(p)A,(q)) 



(1) 



{S(p-q)){G(q)){S(p)) 
where G(q) is the gluon propagator. 



(6) 



In the continuum, the off-shell amputated ver- 
tex function will have the general form 

r M (p 2 ,g 2 ,pg) 

= Fifv + F 2 q fJ , + -FsTm 

+F 4 fa + F 5 i>q^ + F 6 i Vil + F 7 4q^ (7) 
+F 8 ct^P v + F 9 a^q v 

At tree level, this reduces to T° = gojfj,. From 
this we can see that the form factor containing 
the running coupling is F3, while all the other 
form factors are expected to be finite. We can use 
the symmetries of the problem, including current 
conservation, to extract Fs(q 2 ;p = 0). We then 
define the renormalised coupling as 

g R (q 2 ) = Z^Z^Fsiq 2 ) (8) 

From this, we can compute gs.{q 2 ) at various 
(3 values, and relate this to the running coupling 
calculated in other schemes, eg, g p s (q 2 ) or the 
lattice 3-gluon coupling [2]. 

4. Computation and results 

We have analysed propagators generated from 
36 configurations at (3 = 6.0 with the SW ac- 
tion, with a lattice size of 16 3 x 48, at k = 
0.1432,0.1440,0.1450 (corresponding to m PS a = 
0.386, 0.311 and 0.257 respectively). The propa- 
gators have been fixed to landau gauge, with ac- 
curacy yjj- Y^, X \d^,A^(x)\ 2 < 10" 6 , and Fourier 
transformed. 

4.1. Propagator in time 

A complete analysis has been done of the prop- 
agators on 9 configurations. The values for the 
scalar and vector masses are given in table 1. 

Table 1 



Values for the scalar and vector mass of the quark 
at (3 = 6.0, for different values of k. 
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Figure 1. Scalar and vector propagators in time, 
for 36 configurations, at k = 0.1432 

This gives a value for m q in the chiral limit 
of 340 ± 40 MeV, in good agreement with the 
350 ± 40 MeV reported in [1]. For propagators 
with finite spatial momenta inserted, the signal 
deteriorates rapidly with increasing momentum, 
so no meaningful information can be extracted 
with the current statistics. 

4.2. Propagator in momentum space 
Preliminary analysis has been done on 15 con- 
figurations at all 3 k values. We have computed 
and fitted 

Z i ,(p) = -iTr pT(p)/p 2 
and 

m(p) = TiT(p)/Z^(p) 

Zip is shown as a function of pa in fig 2. Z^{jp) has 
been fitted to a logarithmic function Z^{jp) = Z -\- 
k In pa, while ra(y>) has been fitted to a polynomial 
function m(jp) = M + bip + b 2 p 2 , where 62 turns 
out to be consistent with zero. The best fit values 
(setting 62 = 0) are given in table 2. 

4.3. Future work 

The analysis of the quark propagator at (3 = 6.0 
discussed here is currently being extended to in- 
clude study of the vertex function. When this 
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Figure 2. Z^,(jp) as a function of p for 15 config- 
urations at k = 0.1432. The line represents a fit 
to a logarithmic function 



Table 2 

Best fit values for the inverse quark propagator. 



K 


Z 


k 


M 


b 1 /a 


0.1432 


1.97 


-0.67 


0.077 


-0.243 


0.1440 


1.92 


-0.59 


0.058 


-0.235 


0.1445 


1.87 


-0.51 


0.050 


-0.234 



However, more analytical and computational 
work is required in order to fully understand the 
significance of the results reported here. In par- 
ticular, the pole mass seems to contain a number 
of ambiguities, although on the positive side it 
must be noted that it does decrease with increas- 
ing ac, as required. Firstly, the relation between 
the mass extracted from the fall-off of the prop- 
agator in time and the pole mass in momentum 
space is unclear. Secondly, the linear dependence 
of the pole mass on pa is not well understood. 
Hopefully, more light can be cast on both these 
issues by studies at different (3 values. 
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is completed, the propagator and vertex function 
will be computed also at (3 = 6.2 and 6.4 on hy- 
percubic lattices, in order both to get a proper 
understanding of discretisation and finite size er- 
rors and possible power divergences in the quark 
mass, and to obtain a reliable determination of 
the running coupling. 

5. Discussion and conclusions 

It is encouraging that we get such a clean and 
unambigouous signal for the quark propagator 
with a small number of configurations. This gives 
us confidence that we will be able to extract a 
reasonable signal also for the vertex funtion. In 
particular, we see from fig. 2 that the quark field 
renormalisation constant, which enters into equa- 
tion (8) for the running coupling, can be precisely 
determined. 



